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ASSEMBLY  OF  SYSTEMS  HAVING  MAXIMUM  RELIABILITY 

by 

Cyrus  Derman,  Gerald  J.  Lieberman  &  Sheldon  M  Ross 

[0]  SUMMARY 

The  first  problem  considered  in  this  paper  is  concerned  with  the 
assembly  of  independent  components  into  parallel  systems  so  as  to 
maximize  the  expected  number  of  systems  that  perform  satisfactorily,, 
Associated  with  each  component  is  a  probability  of  it  performing 
successfully.  It  is  shown  that  an  optimal  assembly  is  obtained  if  the 
reliability  of  each  assembled  system  can  be  made  equal „  If  such 
equality  is  not  attainable,  then  bounds  are  given  so  that  the  maximum 
expected  number  of  systems  that  perform  satisfactorily  will  lie  within 
these  stated  bounds;  the  bounds  being  a  function  of  an  arbitrarily 
chosen  assembly.  An  improvement  algorithm  is  also  presented 

A  second  problem  treated  is  concerned  with  the  optimal  design  of 
a  system.  Instead  of  assembling  given  units,  there  is  an  opportunity 
to  "control"  their  quality,  i.e,,the  manufacturer  is  able  to  fix  the 
probability,  p,  of  a  unit  performing  successfully  However,  his 
resources  are  limited  so  that  a  constraint  is  Imposed  on  these 
probabilities.  For  (1)  series  systems,  (2)  parallel  systems,  and  (3) 
k  out  of  n  systems,  results  are  obtained  for  finding  the  optimal 
p's  which  maximize  the  reliability  of  a  single  system,  and  which 
maximize  the  expected  number  of  systems  that  perform  satisfactorily  out 
of  a  total  assembly  of  J  systems, 

1 


1 


INTRODUCTION 


In  a  previous  paper  [l]*  che  authors  considered  the  following 
reliability  problem.  A  system  has  n  different  types  of  components 
Associated  with  each  component  is  a  numerical  value  Let 
{am}  (m  =  12,  ;n)  denote  the  set  of  numerical  values  of  the  n 

components  Let  R(a  j  a  .a  )  denote  rhe  probability  that  the 
system  will  perform  satisfactorily,  l  e  ;  R(a\  a^,  .a11)  is  the 

reliability  of  the  system  Now  suppose  a™  ^  a!'  •  a™  are  J 

components  of  type  m  (m  =  1,2,  .n)  Then  J  systems  can  be 

assembled  from  these  components  Let  N  denote  the  number  of  systems 
that  perform  satisfactorily  N  is  a  random  variable  whose  distribution 
will  depend  on  the  way  the  J  svstems  are  assembled  The  results 
obtained  show  that  if  R(a\  a^,  an)  has  the  properties  of  a  joint 
cumulative  distribution  function  then  of  all  different  ways  in  which 
the  J  systems  can  be  assembled,  E(N)  is  maximized  if  these  J 
systems  have  reliability  R(a^  a^.-  fj  =  1,2.  ,J),  l.e  , 

assemble  the  best  of  each  type,  the  next  best  of  each  type.  .  .and 
finally  the  worst  of  each  type 

The  aforementioned  results  <»re  applicable  r.o  series  systems  of 
independent  components,  but  are  not  applicable  to  parallel  systems  of 
independent  components  Section  2  of  this  paper  treats  the  problem 
of  parallel  systems  of  independent  components,  and  shows  that  an 
optimal  asrembly  13  obtained  if  the  reliability  of  each  assembled 
system  can  be  made  equal  Furthermore  if  equal  reliability  for  each 
assembled  system  is  not  possible,  then  bounds  are  obtained  so  that  the 
maximum  expected  number  of  systems  that  perform  satisfactorily  will  lie 


2 


within  these  stated  bounds,  the  bounds  being  a  function  of  an  arbitrarily 

chosen  assembly  An  algorithm  will  be  presented  which  will  produce 

assemblies  thac  result  in  an  improvement  in  these  bounds,  although 

will  not  necessarily  lead  to  the  optimal  assembly  The  results  are  not 

only  applicable  to  the  aforementioned  problem  but  are  also  applicable 

to  the  assembly  or  paraiie^  systems  of  independent  units,  where  the 

units  in  a  system  are  interchangeable  (only  one  type  of  unit  is  in  the 

system)  The  results  are  also  valid  for  the  interchangeable  and  non- 

interchangeable  cases  when  the  number  of  units  of  a  given  tyne  that 

must  appear  in  a  given  system  is  not  fixed  in  advance 

Section  3  is  concerned  with  a  variation  of  the  assembly  problem 

dealing  with  system  design  Suppose  a  single  svstem  is  to  be 

t  h 

constructed  containing  n  units  Attached  to  the  m  unit  is  a 

ch 

positive  value  ,  which  will  denote  the  probability  that  the  m 

unit  will  perform  satisf ac to c  liy  >  and  these  probabilities  are  assumed 

to  be  independent  The  n  units  are  to  be  "manufactured";  and  the 

manufacturer  has  sufticlent  control  of  his  process  that  he  is  able  to 

produce  at  an  aimed  at  "p"  level;  but  his  resources  are  limited  so 

that  there  is  a  constraint  imposed  on  the  values  of  the  p's,  namely 
n 

^  p  =  A  For  (1)  series  svstemst  (2)  parallel  svstems  and  (3)  k 
*-  m 

m=l 

out  of  n  systems  (the  system  operates  satisfactorily  if  at  least  k 
out  of  n  units  operate  satisfactorily),  the  problem  considered  is  to 
find  the  desired  p’s  so  as  to  maximize  the  reliability  of  the  system 
For  a  series  system  the  optimal  solution  is  to  make  aij  the  p’s 
equal  to  A/n  For  the  parallel  system  the  optimaj  solution  is  to  make 
one  of  the  p  s  equal  to  A  and  the  rest  0  For  the  k  out  of  n 


s 


system,  the  results  are  more  complex 


Section  4  is  concerned  with  the  same  problem  considered  In  Section 
3„  except  that  nJ  components  are  to  be  "manufactured"  and  assembled 
into  J  systems  each  being  a  k  out  of  n  system  Again,  it  is 
desired  to  maximize  the  expected  number  of  systems  that  perform 
satisfactorily,  subject  to  constraints  on  the  p  s  A  characterization 
of  the  optimal  solution  is  given 


2„  OPTIMAL  ASSEMBLY  OF  PARALLEL  SYSTEMS  OF  INDEPENDENT  COMPONENTS 

The  following  problem  will  be  considered  first..  A  set  of  M  units 

are  given  which  are  numbered  1.  2.  ,M  The  M  units  are  to  be 

partitioned  into  J  disjoint  systems  After  completion  of  a  partition, 

th 

the  number  of  units  contained  in  the  j  system  (j  -  1,2,  <■ « „  J) 

J  l 

will  be  denoted  by  n  with  the  added  restriction  that  l  n.  -  Me 

J  j*l  J 

A  system  will  perform  satisfactorily  if  at  least  one  of  the  units 

in  the  system  performs  satisfactorily,  i  e  ,  it  is  a  parallel  system 


m 


Attached  to  the  unit,  m  *  1  2  -M.  is  a  positive  value  p 

which  will  denote  the  probability  that  the  mC^  unit  will  perform 
satisfactorily,  and  these  probabilities  are  assumed  to  be  independents 
For  a  given  parition  the  reliability  or  system  j  (j  ■  1, 2  » „  „ , J) , 
j  is  the  probability  that  the  system  will  perform  satisfactorily,  and 


can  be  expressed  as 


1/ 


1 /  A  partition  will  allow  for  one  or  more  systems  to  contain  no  units 
J 

so  long  as  l  n  *  M  The  reliability  of  a  system  containing  no  units 

J-l  J 

will  be  taken  to  be  zero 


rr  i  -it-  -“ruliiifii 


i  -  7T  (i  -  p  ) . 

all  m  in 
system  j 


(1) 


Let  N  denote  the  number  of  systems  that  perform  satisfactorily,  so 
that  N  Is  a  random  variable  whose  distribution  will  depend  upon  a 
given  partition.  For  a  given  partition,  the  expected  number  of 
systems  that  perform  satisfactorily,  E(N),  is  then  seen  to  be 


J  _  J 

e(n)  -  l  [i  -  7T  <i  -  p  )]  -  j  -  l  7T  (i  -  p  >»  (2) 

j»l  all  m  in  j-1  all  m  in 

system  J  system  j 


The  problem  treated  in  this  section  is  to  find  the  partition  which 
maximizes  (2),  or  alternatively,  to  find  the  partition  that  minimizes 
the  expression 


J 


I 

j-1 


7T 

all  m  in 

system  j 


a. 

m 


(3) 


where  a  '■  (1  -  p  ),  Henceforth,  the  a  will  be  referred  to  as  the 
mm  m 

positive  value  attached  to  the  m  unit. 

In  a  given  partition,  denote  the  set  of  n^  units  appearing  in 
system  J  (j  ■  1,2,,,„,J)  by  A^  and  define  | A ^  |  to  equal  the 
product  of  the  values  of  the  units  in  the  set  A^.  If  A^  is  the  null 
set,  |Aj|  is  defined  to  be  1,  Define 

J 

T  -  Minimum  £  jA.  |  ,  (4) 

J-1  J 


5 


vh.r.  th.  .into  1.  ton  ovnr  .11  fon.ibl.  partition,  of  tho  M 
unit,  into  J  dlajolnt  ..«•  Th.  major  ...ult  of  thl.  ..ctlon  .how. 
that  for  any  partition  ^  A,. . . .  ^ .  t  ..tl.fl..  th.  following 

inequalities, 


j 

J  Min  |A  |  1  T  <  l  |A, 

1  j-1  3 


In  order  to  prova 


this  result,  two  lemmas  must  be  verified. 


Lemma  1»  Let  a 


aud  b  (m  -  1,2 . s>  be  poaitive  numbera, 


a  ® 

7f  a  -  7T  bm, 

m  m  '*  m 

m-1  m-1 


l  (a/b)  i  s. 
*..mm 


Proof  l  Sint.  .  /b.  -  TT  b  /  TT  *„.  th.  contlo.lon  of  th.  1«— 
88  m-1  m-1 

will  follow  if  it  can  be  shown  that 


3-1  S-l 


l  +  TT  b  /  7T  am  i  8‘ 

n.l  m  °  m-1  m  m-1 


Now,  define  the  function 


s-l  s-l 


£<Vi . Vi>  ■  i,  ( VV  +  £  V  &  *. 


Differentiating  this  function  with  respect  to  b^,  yields 


Sf(b1,b2,...,ba-1) 

3  b. 


„  s-1  s-1 

-a. /b.  +  Tf  b  /  Tf  a  ,  k  ■  1(2( , . ,  ,S"1  ■ 
K  K  ,  m  -  O' 

m«l  m*l 


Equating  these  to  zero,  results  in 


-  b  7  rr 


an/bi,  -  77  b  /  77  a 

k  k  .  m  ,  m 
m**l  ra»l 


a  -  c,  k  -  1,2, ....s-l. 


This  set  of  equations  implies  that  each  a^/b^  equals  the  constant,  c. 
If  all  of  the  (s-1)  terms  of  the  left  hand  side  of  (8)  are  multiplied 
together,  then  it  follows  from  (8)  that 


7T  <Vbk>  *  c 

k-1 


(s-1) 


However,  from  (8),  the  left  hand  side  of  (9)  is  also  equal  to  c.  Hence 
c  must  equal  1. 


Since  f(b^,  b2,..„,bg  approaches  infinity  as  one  or  more  of  the  b's 

approach  zero  or  infinity,  f  (b^  ,b2, .  „ ,  ,bg_^)  is  minimized  when  its 

partial  derivatives  vanish.  Thus,  it  attains  it  minimum  when  a  ■  b 

mm 

(m  -  1,2,... , s-1),  and  f  (b^,b2,. »  .  *bg<_^)  >_  s  so  that  the  lemma  is 
proved. 


Lemma  2s  Let  A^,  A2,„„u,Aj  be  a  partition  such  that 


I A1  i  “  1^1 

optimal. 


|Aj  j .  Then  T 


i.e.,  this  partition  is 


Proof:  Consider  any  other  partition  -  say  B1,  B 
theoretic  identities,  it  follows  that 


Using  set 


|Aj  bj 


where 


Aj  and  B^  are  the  complements  of  and  B  ^ , 
to  prove  the  lemma  it  is  sufficient  to  show  that 


respectively.  Hence, 


Inequality  (10)  follows  from  Lemma  1  by  noting  that 


7T 

J-l 


iBJAJ 


J 

Tf  |A  B  I, 
j*l  3  J 


because  both  terms  are  equal  to 


(10) 


i(U  AB)j, 
J-l  3  3 


Theorem  1. 

For  any  partition  A£ 


,AJt  T  satisfies  the  following 
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inequalities 


J  Min  I A 


^li'i  I  ia,i- 

J  j  =  l  3 


Proof : 


Suppose  that  |A^|  »  min|Ajj.  It  then  follows  that  the  units  in 


the  sets  A.,  1  +  k  have  values,  a  ,  associated  with  them  such  that 

j  m 


| Aj  |  l\l*  Nom»  consider  a  new  set  of  M  units.  These  units  are 

partitioned  in  the  same  manner  as  the  original  units.  The  values 
associated  with  the  units  in  set  A^  are  the  same  as  those  in  the 
original  set.  However,  the  values  associated  with  the  units  In  set 


Aj,  j  k  are  all  less  than  or  equal  to  the  corresponding  a's  in 


the  original  set,  but  are  such  that  j A^  (  now  equals  jAjJ  for  all 


j  y  k.  For  this  new  set  of  M  units,  let  T_  denote  the  minimum 
J 


\  |A  |  ,  where  again  the  minimum  is  taken  over  all  possible  partitions 

J-l  J 

of  the  M  units  into  J  disjoint  sets.  However,  from  Lemma  2,  the 


given  partition  must  be  the  optimal  one  for  the  new  set  of  units,  and 
furthermore,  Tq  *  J" |  j .  Since,  the  function  T  is  obviously  a 
monotone  increasing  function  of  the  set  of  values  {a^,  it 

follows  that 


T  >  Tq  =  jJaJ  =  J  minjAjj. 


The  other  inequality,  i.e.,  T  _<  \  j  A .  |  ,  is  obvious.  Hence,  the 

j-l  J 


theorem  is  proved. 

The  theorem  indicates  that  the  maximum  expected  number  of  systems 
that  perform  satisfactorily  will  lie  within  the  stated  bounds,  these 

9 


. .  r  -r  - n 


bounds  being  a  function  of  the  chosen  partition  Furthermore,  if  a 
partition  can  be  found  that  makes  each  system  have  the  property  that 
the  product  of  the  probabilities  of  each  unit  failing,  i.e.,  the  a's, 
equal,  then  this  partition  is  optimal  in  that  E(N)  i9  maximized. 

These  results  lead  to  questioning  whether  or  not  an  algorithm  can  be 
obtained  which  will  determine  the  optimal  partition.  Unfortunately, 
the  authors  have  been  unable  to  find  one,  but  an  algorithm  will  be 
presented  which  should  lead  to  a  "good"  solution;  but  not  necessarily 
an  optimal  one. 

A  given  partition  results  in  a  sequence  of  sets  A^,  A2»...,Aj 
(each  set  representing  a  system)  and  a  corresponding  j A^|  ,  j  A^  j  » « • . ,  |  A j 
It  can  be  assumed  that  the  j  A |  * s  are  not  equal;  otherwise,  an  optimal 
partition  has  been  obtained  Choose  any  two  systems  whose  j  A |  's  are 
not  equal  and  without  loss  of  generality,  denote  them  by  A^  and  Aj, 
with  | A2 j  >  jA^j.  It  will  be  shown  that  under  certain  conditions 
units  of  one  can  be  interchanged  with  units  of  the  other,  thereby 
resulting  in  a  new  partition  with  sharper  bounds  than  given  In  (5).  Let 
denote  the  product  of  the  values  which  are  attached  to  those  units 
which  are  to  be  removed  from  A^  and  placed  Into  A2.  Similarly,  let 
«2  denote  the  product  of  the  values  which  are  attached  to  those  units 
which  are  to  be  removed  from  and  placed  into  A^.  Thus,  in  the 

new  partition,  the  |A|’s  are  given  by  (o^/oi^)  j  [  ,  (u^/'o^)  |  ^  j  , 

|  A3 1 . I  Aj  j  -  The  main  result  to  be  obtained  is  that  if 

i  A2  i  >  jAjJ  and  -  (:y  \\\ |  -  |  i  A2  i  “  I  Ai  1  fi  then 

1/  The  symbol  |jz|  reads  the  absolute  value  of  Z, 
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|Aj  +  j  A2  j  >  (a2/ai)|A1|  +  (a1/'u2)|A2i. 


(ID 


Before  proving  this  result,  two  lemmas  are  required 

Lemma  3:  If  |  A2  j  >  (aJ,  then  |  A1  j  +  j  A2  j  '  j  Aj  +  (u^/a^jA^ 

<!*>  |A1j/|A2j  <  '  1. 

Proof:  It  is  evident  that  |  A^  j  +  j  A2  j  >  |  A.  j  +  (u./ylAj 

<">(a2  -  a^) [ j A2 j /a2  -  j A^ j /a^j  >  0.  If  jAj/jA^  -  <  1,  then 

(a2  -  “i)  1 1 Ag  j  /a2  -  iAj/aJ  >  0, 

Now  suppose  that  (a2  -  a^)[|A2j/a2  -  |A^|/o^]  >  0.  This  implies  that 
if  a2  -  <  0,  then  l|A2|/'a2  -  ] A^ j  /a 0,  However,  this  cannot 

hold  because  a2  -  0  and  jA2|  '  |A^|  implies  that 

[  | A2  j  /c»2  -  |A^j/a^]  -  0.  Hence,  (a2  -  <y  '  0  and  consequently, 
[|A2j/a2  -  jAj/a^]  '■  0.  Therefore,  j A^ |  / j A2  j  «•  a^/a2  *-  1,  and  the 
Lemma  is  proved. 

Lemma  4:  If  | A2 J  jAj  and  j Aj^ j  +  { A2 j  '  (oy'yjAjj  +  (otj/o^)  j A2|  , 

then  |  |  A2  j  -  j  Ax  1 1  '  |  (a±/a2>  j  A2  j  -  (o^/yjA-JI 

Proof:  Define  the  function  F(x)  ■  x|a2)  -  (l/x)jAjj,  0  £x£  1, 

end  note  that  it  is  monotone  increasing  from  F(0)  -  -  »  to 

F(l)  *  |  A2 1  -  j  Aj^ ( ,  with  F(|A1j/jA2j)  -  jAj  -  |a2|.  Therefore,  for 
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all  x  such  that  IaJ/IAjI  <  x  <  1, 


|x|A2I  -  (1/X)|A1I|  <IiA2l  -  <12) 

However,  from  Lemm.  3,  ^i/i^l  V«2  <  1.  Thu.,  inequellty  (12) 
la  satisfied  for  x  -  uj/«2.  a”d  the  lemma  is  proved. 

The  converse  of  Lemma  4  Is  stated  as  Theorem  2. 

Theorem  2: 

If  |  A2  i  -  |AjJ  and  | (ai/a2> I A2 I  ”  (ot2 /^i>  1  Ai !  |  iAllf 

then 

|ax|  +  1  a2  i  '  (a2/a1>lA1i  +  ('i1/°2)lA2l  ’ 

Proof:  From  the  monotonicity  of  the  function,  F,  defined  in  the 

proof  of  Lemma  4,  and  the  values  of  F(  j  Ax  j  /  i  A2 1 )  and  F(l)  given  in 
the  proof  of  Lemma  4,  it  follows  that  (12)  is  satisfied  for  only  those 

values  of  x  such  that  IaJ/^I  <  *  <  1-  i8  then  clear  that 

| A  j/|A  I  <  0l/a2  <  1  is  implied  by  the  hypothesis  of  the  theorem. 

The  conclusion  of  the  theorem  follows  from  Lemma  3. 

Corollary  1: 

For  a  given  partition,  if  |A ,j  is  the  minjA^,  J  -  1.2 . J» 

and  lAjj  is  any  other,  1  Aj j  -  i Aj_ j .  and  V  ^  are  such  that 
Theorem  2  holds,  then  mini (o^/c^)  |  A^  j  ,  a2 )  I  Aj  i  ’>  Ai i  This 


12 


implies  that  the  new  partition  results  in  a  higher  lower  bound  than 
that  given  in  (5) . 

Proof : 

From  the  conclusion  of  Theorem  2,  jA^J  +  j  j  (c^/a^lAjJ  +  (c^/o^)  |A2 
But  from  Lemma  3,  this  Implies  that 

i  |  / 1 A2  |  »•  a1/a2  <  1. 

Now,  min{  (o^/a^)  |  A^j  ,  (a^/a^)  j  Aj  |  }  *  j  A^ |  if  and  only  if 

^a2^al^^J  '  I ^i  I  an<*  (o-i/ <^2 ^  I ^2  i  '  1^1  i*  But  latter  require¬ 

ments  are  precisely  the  inequalities  of  Lemma  3,  i.e.  ,  j A^ | /  j A2 1  '-c^/c^  < 

Theorem  2  and  Corollary  1  show  that  one  iteration  can  sharpen  both 
the  upper  and  lower  bounds  given  by  (5).  In  fact,  heuristically ,  one 
seeks  partitions  that  tend  to  equalize  the  |A|  s,  and  this  can  be  done 
systematically  by  interchanging  units  from  one  system  with  units  from 
another  system;  these  units  satisfying  the  conditions  ol  Theorem  2, 
e.g.,  interchanging  units  within  the  highest  and  lowest  j A| * s .  The 
algorithm  would  be  continued  until  there  are  no  pairwise  interchanges 
satisfying  the  conditions  of  Theorem  2.  When  this  occurs,  the  solution 
is  "good"  but  not  necessarily  optimal.  This  can  be  seen  by  examining 
the  following  counterexample.  Suppose  there  are  9  units  with  associated 
a‘s  to  be  divided  into  three  systems  as  follows: 
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.00013527 


System  1: 

.09, 

.0501, 

.03; 

1  Ax|  - 

.00013527 

System  2: 

.06, 

.20, 

.02; 

i  Aj  1 

.000240 

System  3: 

.12, 

.1001, 

.015; 

1  A3  j  - 

.00018018 

For  this  partition 

iAxi 

+  |a2|  + 

|a^[  ■  .00055545,  and  cannot 

be  Improved  by  pairwise  interchanges 

However, 

consider  the  following 

partition: 

System  1: 

.12, 

.0501, 

.03; 

IaJ  - 

.00018036 

System  2: 

.06, 

.20, 

.015; 

|Aj|  - 

.000180 

System  3: 

.09, 

.1001, 

.02; 

i  a3  j  - 

.00018018 

For  this  partition. 

which 

is, 

in  fact, 

optimal, 

iAj  +  |a2| 

+  |A3|  -  .00054054. 

Since  the  aforementioned  algorithm  only  allows  for  pairwise  interchanges, 
it  need  not  lead  to  optimal  solutions . 

The  problem  considered  to  date  has  been  in  the  context  of  taking  M 
units  and  partitioning  them  into  J  disjoint  systems,  with  units  being 
Interchangeable  and  the  number  required  for  each  system  not  specified. 
Suppose  the  problem  is  now  changed  so  that  each  system  contains  n 
different  (non-interchangeable)  types  of  components  and  J  systems  have 
to  be  assembled  from  the  M  •  nJ  units,  i.e.  there  are  J  units  of  each 
type  available.  Again,  each  system  operates  as  a  parallel  system,  and 
as  before  the  objective  is  to  find  the  partition  that  maximizes  the 
expected  number  of  systems  that  perform  satisfactorily.  How  does  this 
new  problem  compare  with  the  problem  previously  treated?  Fortunately, 
the  admissible  partitions  for  this  new  problem  is  a  subset  of  the 


partitions  of  the  original  problem,  and  furthermore,  none  of  the  results 
depended  on  the  number  of  units  assigned  to  each  system.  Hence,  all 
the  results  previously  obtained  are  applicable  to  the  new  problem.  Of 
course,  similar  comments  can  be  made  for  the  case  of  interchangeable 
components,  but  where  each  system  must  contain  n  =  M/J  units. 

3.  SINGLE  SYSTEM  PROBLEM 

Another  variation  of  the  assembly  problem  is  concerned  with  system 

design.  Suppose  a  single  system  is  to  be  constructed  containing  n 

units.  Three  cases  will  be  considered,  namely  (.1)  the  system  will 

perform  satisfactorily  if  all  of  the  n  units  performs  satisfactorily, 

i.e. ,  It  is  a  series  system;  (2)  the  system  wiil  perform  satisfactorily 

if  at  least  one  of  the  n  units  performs  satisfactorily,  i.e.,  it  is 

a  parallel  system;  and  (3;  the  system  will  perform  satisfactorily  if 

at  least  k  >  1  units  periorms  satisfactorily,  i.e  ,  it  is  a  k  out 

ch 

of  n  system.  Attached  to  the  m  unit  is  a  value  p^,  m  -  l,2,...,n, 

which  will  denote  the  probability  that  the  mth  unit  will  perform 

satisfactorily,  0  ^  i,  and  these  probabilities  are  assumed  to 

be  independent.  The  n  units  are  to  be  "manufactured",  and  the 

manufacturer  has  sufficient  control  of  his  process  that  he  is  able  to 

produce  at  an  aimed  at  "p"  level,  but  his  resources  are  limited  so 

that  there  is  a  constraint  imposed  cn  the  values  of  the  p's,  namely 
n 

\  p  =  A,  where  A  is  a  fixed  positive  number.  The  problem  is  to 
m=l 

find  the  desired  p's  so  as  to  maximize  the  reliability  of  the  system. 
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Case  1:  Series  System 

For  a  series  system,  the  problem  is  to  find  the  p^,  P2»,,,»Pn 


which  maximizes  the  reliability,  R,  where 


n 

IT  P_  » 

i  m 

m«l 


subject  to 


0  <  p  <1,  m  ■  l,2,...,n  and  Y  p  ■  A. 
—  m  —  L,  *m 

m-x 


It  can  be  assumed  that  A  <_  n,  otherwise  the  solution  is  to  choose 

n 

each  p^  -  1.  The  problem  is  equivalent  to  maximizing  ][  log  p^ 


m-1 


subject  to  the  same  conditions.  Since  log  p^  is  a  concave  function 
and  ignoring  the  constraints  Pm  £  1,  m  »  l,2,...,n,  it  is  well  known 
that  the  optimal  values  of  p^,  P2»”,*Pn  are  Pj  “  P$  P{{  “  A/n. 

Since  A  ^  n,  the  ignored  constraints  are  satisfied.  Hence 
pf  *  p^  * • • •■  p*  =  A/n  is  optimal. 


Case  2:  Parallel  System 

For  a  parallel  system,  the  problem  is  to  find  the  p^,  P2»,,,*Pn 
which  maximizes  the  reliability  R,  where 


r  .  i  -  7T  (i  -  pn>, 

m*l 


subject  to 


0  <_  p^  _<  1,  m  ■  1,2,. . .  t,n  and  £  pm  *  A. 


m«l 


m 
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It  can  be  a _ -ed  that  a  i,  otherwise  the  optimal  solution  is  to 

choose  at  least  one  ot  the  p's  equal  to  1  and  the  rest  arbitrary  (but 

n 

subject  to  the  constraint  that  £  p  »  A) .  The  problem  is  equivalent 

m*i 

to  minimizing 


n 

1  Aog(l  -  pm), 


subject  to 

n 

pm  >  0,  m  *  l,2,  -.,n,  and  V  -  A. 

m»i 


The  constraint  p  =1,  m  *  1,2,...  ,n  Is  superlluous  since  A  '  1. 

m  n 

The  functions  log(l  -  p  )  are  concave.  Therefore  \  log(l  -  p  ) 

m“i 

is  a  concave  function  m  (p^,  P2»*  -»Pn)  with  the  minimum  at  an 
n 

extreme  point  of  £  p  *  A,  p  ^  0  (m  ■  1,2,.  ,n).  Each  extreme 
m-1 

point  has  one  p  ■  A  and  the  remaining  equal  to  zero  Thus,  in 
m 

particular,  pj  *  A,  p£  =  p^  »  « p£  «  0  is  optimal. 


Case  3:  k  out  of  n  System 

For  a  k  out  of  n  system,  the  system  wili  perform  satisfactorily 

if  at  least  k  out  ot  n  units  perform  sauistacto-ily .  Let  X  , 

m 

m  *  l,2,..,n,  be  independent  Bernoulli  random  variables  with  parameter 

p  ,  Let  Y  ■  X,  +  X„  X  .  Then  the  reliability  R  can  be 

m  1  z  n 

expressed  as 


R  =  P{Y  >  k}. 
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The  problem  Is  to  maximize  R  subject  to 

n 

'i'.-1'  and  l  pm  -  *• 

m-1 

It  can  be  assumed  that  A  <  k,  otherwise  the  optimal  solution  Is  to 
choose  k  of  the  p's  equal  to  1  and  the  rest  arbitrary.  If  the 
constraints  pm  <_  1  are  ignored,  the  Kuhn-Tucker  conditions  indicate 
that  if  pj,  p^.-.-.p*  are  optimal,  then  there  exists  a  number 
X*  _>  0  satisfying  the  following: 

If  p*  *  0,  then  — —  PtY  s  k}  -  A*  <  0  at  p  *  p*, 

rm  dp  —  —  m  m 

*ra 

for  m  *  1,2,. . ,  ,n. 

If  p*  >  0,  then  -r^~  P{Y  >  k}  -  X*  »  0  at  p  -  p*, 

m  dp  —  mm 

rm 

for  m  ■  1,2,, .  ,n, 

n 

If  A*  '  0,  then  T  p*  -  A- 

m»l 

Suppose  there  are  r  variables  which  are  greater  than  zero  in  the . 
optimal  solution  and  denote  them  by  pj,  pg,...,p*.  The  P{Y  _>  k}  can 
be  expressed  as 

P{Y  >_  k}  -  p^^  F(p2,  p3, .  ,pn>  +  U  -  p^)  G(p2,  P3, .  •  •  »Pn) » 

where  the  function  F  is  P{X0  +  X-  +••  +  X  >  k  -  1}  and  the  function 

Li  n  — 


G  is  P{X2  +  X3  +•••  +  Xn  >  k}. 


Hence , 


•— —  P{Y  k}  *  F(p2 ,  P^»’«*»Pn)  —  G(p2 ,  p3»..»»Pn)  *  H(p2 ,  P3»*>*>Pn)> 


Therefore,  H(pi|,  p*,...,p£)  -  A*  can  be  solved  uniquely  for  p*  as 
a  function  of  p^,  p*,,».,p*  and  A*, 

Similarly, 


—  P{Y  >  kl  •  H(Pl,  p3 . pn) 


can  be  solved  uniquely  for  p*  as  the  same  function  of  p*,p£ . p* 

and  A*  so  that  pj  =  p|.  In  the  same  manner,  it  can  be  shown  that 

p*  «  p*  p*  ,  and  furthermore,  it  is  clear  that  p*  *  A/r, 

i  *  l,2,.,.,r.  If  ^  <  1,  then  the  ignored  constraints  are  satisfied. 

This  is  easily  shown  to  be  the  case.  It  is  clear  that  since  if 

more  than  (n  -  k)  p’s  are  zero,  then  PiY  >  k}  -  0,  which  cannot 

A  A 

be  a  maximum.  Therefore  -  £  £  1,  so  that  the  ignored  constraints 

are  satisfied.  The  foregoing  results  do  not  indicate  how  to  determine 
r.  One  way  is  to  evaluate  P{Y  >_  k}  with  p*  *  A/r,  i  *  l,2,...,r, 
choosing  that  value  of  r  »  r*  which  maximizes 


•  £  WHI 


f ][— F  (l-  — lr~^  for  r»k,  k+l,,..,n.  (13) 
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Some  insight  is  obtained  by  replacing  the  right  hand  side  of  (13) 
by  the  normal  or  Poisson  approximation,  depending  upon  which  is 
appropriate  When  r  is  large  and  A/r  is  near  0  or  1  the  Poisson 
is  appropriate;  when  A/r  is  "near"  1/2,  the  normal  is  appropriate. 
For  either  case,  k  (and  consequently  r)  and  n  should  be  large. 

For  the  normal  approximation  r  a  r*  is  to  be  chosen  which  maximizes 


✓r(A/ r) (1-A/r)  /A(l-A/r) 


For  r  =  k,  k+l,,..,n  Since  k-A  ^  0,  clearly  r*  »  n.  Thus, 

P*  “  ?2  Pn  “  A//n'  (15) 

If  the  Poisson  approximation  is  used  for  the  smaller  values  of  A,  the 
right  hand  side  of  (13)  yields 


P{Y  k)  =  1  -  l  e~AAJ  , 

j=0  J  ‘ 

which  is  Independent  of  r.  That  is,  approximately  speaking,  the 
choice  of  r  has  very  little  influence  on  the  left  hand  side  of  (13). 
Therefore,  it  can  be  said  that  if  n  and  k  are  both  large,  an 
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approximately  optimal  solution  is  given  by  (15) 


The  results  given  in  (15)  are  consistent  with  the  n  of  n  case 
(series  system).  It  does  not  appear  to  be  consistent  with  the  1  of 
n  case  (parallel  system).  However,  this  is  not  unexpected  since  k 
must  be  large  in  order  tor  the  approximation  to  be  good. 

The  solution  given  in  (15)  is  easily  shown  to  be  exact  and  optimal 
if  the  problem  considered  is  tc  maximize  the  variance  of  the  random 
variable  Y,  the  number  of  units  that  perrcrm  satisfactorily ,  subject 
tc  the  usual  constraints  (Since  A  k  and  A  is  the  expected 
number  ot  components  that  function  satisfactorily,  maximization  of  the 
variance  is  desired  ) 

4.  MULTIPLE  SYSTEM  PROBLEM 

The  problem  considered  in  this  section  is  the  same  as  that 
considered  in  Section  3,  except  that  nJ  components  are  to  be 
"manufactured"  and  assembled  into  J  systems.  Each  system  performs 
satisfactorily  if  at  least  k  out  or  n  units  perform  satisfactorily. 
It  is  desired  to  maximize  the  expected  number  ot  systems  that  perform 
satisfactorily,  E(N),  subject  to  constraints  on  the  p's 

Motivated  by  the  results  obtained  in  Section  3  using  the 
approximation,  It  will  be  assumed  initially  that  the  probability  of 
each  unit  performing  satisfactorily  within  the  system  will  be  equal, 
i.e.,ali  units  in  the  jC^  system  have  probability  p^  of 
functioning  The  problem  is  to  determine  p^,  p^,.,  ,Pj  to  maximize 

D(Pi  ,  p2,...,pj  -  ll  l^lpfu  -  P^n_1  (16> 

j-i  i=k'  J  J 
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subject  to 


OlPjil.  J  -  1.2 . J 


I  1  A» 

j-1  J 


where  0  <  A  <  nJ. 


The  Kuhn-Tucker  conditions  indicate  that  if  p*,  p*,...,p*  are 
optimal,  then  there  must  exist  numbers  X*  _>  0,  u*  j>  0,  j  -  l,2,...,n, 


satisfying  the  following: 


if  pj 


g 

0,  then  D(pltp2, . . .  .pjJ-nX*  -  y*  _<  0  at  Pj 


for  J  m  1,2,. . . ,J. 


(17a) 


If  p*  >  0,  then  D(p1,p2,...  ,pj)-nX*  -  y*  »  0  at  p^  - 


for  j  -  1,2,. .. ,J. 


(17b) 


If  v*  -  0,  then  P*<1 


for  j  ■  1,2, ... ,J. 


(17c) 


If  y|  >  0,  then  p*  -  1 


for  j  -  1,2 . J- 


(17d) 


If  X*  -  0 


si 

,  then  n£  p?  <_  A. 

j-1  J 


(17e) 
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C17f) 


If  X*  >  0, 


Chen 


J 

nI  P? 
j-1  J 


A  • 


Now, 


3D 

3pj 


n-li  k-1 

n  n-ki  Pj 


(1 


V 


n-k 


3D 

so  that  - —  >  0  for  0  <  p  <1  and  equal  to  zero  when  p  ■  0  or  1. 

J  —  J  J 

J 

If  p*  *  0«  then  p*  a  1,  When  p$  ■  1,  - —  *  0  sc  that  v*  ■  -nX* 

J  J  J  J 

from  (17b);  this  is  a  contradiction,  and  hence,  v*  ■  0  for  all  j, 

j  •  1,2 . J.  Now  if  any  p*  *  1,  then  from  (17b)  X*  »  0,  so 

that  every  p*  must  be  zero  or  one  (from  17a  or  17b)  since  ■  0 


only  at  these  values  of  p.  If  0  <  <  1,  then  from  (17b) 


j 


n(n-k]  pj(1  "  Pj)n_k  "  nX*  >  °‘ 


For  a  given  X*  the  equation  has  at  most  two  possible  roots 
range  0  <  p^  <  1.  Therefore,  an  optimal  solution  is  of  the 

pj  -  p$ - P$  ■  o,  and  pjrt  -  p;t2 - P,rty  -  p. 

Px+y+l*  Pjc+y+2 » •  ■  •  »p j  ■  p»  for  80me  x  and  y»  where  P*  P 


in  the 
form 

are 


solutions  to 


n  j"  (1  -  p)n  ^  ■  X*,  for  some  X*  >  0  (18a) 

and  vp  +  (J  -  x  -  y)p  -  ^  .  (18b) 

An  optimal  solution  may  consist  of  [s]  p’s  each  having  value 
one  and  (J  -  [s])  p's  each  having  value  zero,  where  [s]  is  the 
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greatest  Integer  less  than  or  equal  tc  A/n  But  this  is  consistent 


with  the  given  form  when  y  =■  0,  p  *=•  0,  p  ■  1, 
n  ^  ^ 

The  curve  £  [  ” Jp  (1  -  p)n  being  s  shaped  is  concave  over 

Is  k 

k-1 

the  range  — r-  ^  p  <_  1,  It  is  evident  that  if  A/nJ  satisfies 

k™  X  A  A 

— -  <  — -  '  1,  then  p*  *  p*  1 “  p1?  *  — r  is  optimal, 

n-l  nJ  i  *  nj 

When  k  «  n,  the  series  case,  the  optimal  solution  is  to  have  p^l 
for  j  «  1,2,.  ,[A/nJ,  P^A/n]+1  *  A/n  -  [A/n],  and  p^  *  0  for 
j  *  [A/n]+2,  ,J  This  is  immediate  from  convexity  arguments, 
independent  of  the  foregoing  argument. 

At  the  outset  of  this  section  it  was  assumed  that  the  probability 
of  each  unit  performing  satisfactorily  within  a  system  was  equal.  This 
may  not  always  be  reasonable  In  Section  3,  it  was  evident  that  when 
A  ^  k ,  the  optimal  solution  was  to  choose  k  of  the  p's  equal  to  1 
and  the  rest  0  In  the  problem  encountered  in  this  section  where  nJ 
units  are  to  be  assembled  into  J  systems,  it  may  very  well  be  that 
allocating  k  of  the  tocal  resources  to  a  system  is  appropriate  even 
if  A  _  kJ  and  hence,  the  initial  assumption  about  equal  probabilities 
within  a  system  may  not  be  valid  Thus,  any  algorithm  which  seeks  an 
optimal  solution  tc  this  problem  requires  a  consideration  of  such 
possible  systems,  and  their  contribution  to  the  total  of  the  expected 
number  of  systems  that  perform  satisfactorily  as  well  as  allocating  the 
remaining  resources  so  as  to  maximize  D  in  equation  (16). 


REFERENCES 


Deraan,  C. ,  Lieberman,  G.J. ,  and  Ross,  S.M.,  On  Optimal  Assembly 
of  Systems,  Technical  Report  No.  140,  Department  of  Operations 
Research  &  Department  of  Statistics,  Stanford  University, 

August  25,  1971.  (To  appear  in  Logistics  Research  Quarterly 
December  1972.) 


